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The Kondo effect in a one-dimensional spin-i XXZ model in the gapless XY regime ( — 1<A<1) 
is studied both analytically and numerically. In our model an impurity spin (5 = 1/2) is coupled to 
a single spin in the XXZ spin chain. Perturbative renormalization-group (RG) analysis is performed 
for various limiting cases to deduce low-energy fixed points. It is shown that in the ground state 
the impurity spin is screened by forming a singlet with a spin in the host XXZ chain. In the 
antiferromagnetic side (0 < A < 1) the host chain is cut into two semi-infinite chains by the singlet. 
In the ferromagnetic side ( — 1 < A < 0), on the other hand, the host XXZ chain remains as a single 
chain through "healing" of a weakened bond in the low-energy (long-distance) limit. The density 
matrix renormalization group method is used to study the size scaling of finite-size energy gaps and 
the power-law decay of correlation functions in the ground state. The numerical results are in good 
agreement with the predictions of the RG analysis. Low-temperature behaviors of specific heat and 
susceptibility are also discussed. 
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I. INTRODUCTION 

There has been recent resurgence of interest in 
the Kondo efFectS in one-dimensional (ID) strongly- 
correlated systems. In ID interacting systems belong- 
ing to the universality class of the Tomonaga-Luttinger 
(TL) liquids, a static impurity potential has a drastic ef- 
fect and is renormalized to infinity or zero, depending on 
whether the interaction is repulsive or attractiveE This 
anomalous response to a static impurity of TL liquids 
has attracted a lot of attention and led to further studies 
on effects of a dynamic impurity (typically a magnetic 
impurity) in a TL liquid. A generalized Hubbard model 
with an impurity spin {S = ^) and its variants have 
been studied by many authors. It was found that the 
Kondo temperature, which is a typical energy scale for 
host electrons to screen an impurity spin, has a ByOjwer-law 
dependence on the Kondo exchange coupling.Bcl Prop- 
erties of low-energy fixed points have been discussed us- 
ing perturbative renormalization group analMsisu and the 
boundary conformal field theory approach.Em A recent 
Monte Carlo studya on the susceptibility of an impurity 
spin is consistent with anomalous uower-law tempera- 
ture dependence conjectured earlier .□ In addition to the 
models with a simple Kondo coupling, there are some ex- 
actly solvable models in which an impurity spin is coupled 
to the spin density of electewe via special forms of the 
Kondo exchange couplingJlJ'til The results obtained for 
these models using the Bethe-ansatz technique, however, 
do not completely agree with the previous studies,Q'la and 
this remains as a question to be resolved. 

In this paper we consider a simplified model which 
we believe shares common features with the above- 
mentioned Kondo effect in ID interacting electronic mod- 



els like the Hubbard model. We here focus on the spin 
sector and discard the charge degree of freedom. This 
may correspond to the half-filled case in the original elec- 
tronic models. The Hamiltonian of the system we discuss 
in this paper has the form H — Hq + Hk, where Hq de- 
scribes the host S ~ ^ XXZ spin chain, 

Ho = jY^ {SfSf^, + S^Sy^, + AS^S^,) , (1) 

i 

and Hk the Kondo coupling, 

The size of the impurity spin is also assumed to be 1/2. 
An important point of our model is absence of SU(2) 
spin rotation symmetry. We assume |A| < 1 to ensure 
that the host XXZ spin chain has gapless excitations. 
For simplicity we have used the same parameter A in 
Hq and Hk- We note that A in Hq is an important 
parameter controlling power-law behavior of various cor- 
relations while A in Hk does not play any significant role 
in the following discussions. The Kondo coupling Jk can 
be either antiferromagnetic or ferromagnetic, but we will 
concentrate on the antiferromagnetic case {Jk > 0) in 
this paper. r-irn 

Eggert and AfFleckll3't2l studied, among various kinds of 
disorder, the model at the isotropic point (A = 1). They 
concluded that the impurity spin Simp forms a singlet 
with So, and that the Heisenberg chain is decoupled into 
two semi-infinite chains in the low-energy limit. A lead- 
ing irrelevant operator at the fixed point was identified 
and shown to have scaling dimension 2. It corresponds 
to exchange coupling between boundary spins of the two 
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decoupled chains. In this paper we extend their analysis 
to the XXZ case (|A| < 1). We first bosonize the Hamil- 
tonian and study its renormalization group (RG) flows in 
the weak-coupling limit and in the strong-coupling limit. 
We will argue that the system is renormalizcd to stable 
low-energy fixed points where the impurity spin (5 = ^) 
is screened exactly. At the fixed points the boundary 
condition for the host XXZ spin chain depends on the 
parameter A of the host chain: For < A < 1 the 
spin chain is cut into two semi-infinite chains with open 
boundary condition at i — ±1. On the other hand, for 
— 1 < A < the host spin chain is not affected much 
by the singlet and stays as a single chain. Leading ir- 
relevant operators at these fixed points have noninteger 
scaling dimensions, yielding noninteger power-law tem- 
perature dependence of impurity contribution to specific 
heat and susceptibility. As evidences for this picture we 
will show finite-size scaling of energy gap and spin-spin 
correlation functions in the ground state, both of which 
are obtained by using the density matrix renormaliza- 
tion group (DMRG) method. The numerical results are 
consistent with the picture drawn from the perturbative 
RG analysis. We note that outj-esults are very differ- 
ent from a recent paper by Liujij who studied the same 
model as ours using mysterious transformations and cal- 
culated various quantities near a strong-coupling fixed 
point. For example, he obtained superlinear temperature 
dependence (T": a > 1) for the impurity contribution to 
the specific heat and vanishing susceptibility at zero tem- 
perature, both of which cannot be correct from general 
grounds. 

The plan of this paper is as follows. In section II we 
discuss RG flows of our model using the standard abelian 
bosonization method. Impurity contributions to specific 
heat and susceptibility are also discussed. We show re- 
sults of numerical DMRG calculations in section III and 
compare them with conclusions of the perturbative RG 
in section II. For simplicity we set J = 1 throughout this 
paper. 



II. PERTURBATIVE RENORMALIZATION 
GROUP ANALYSIS 

A. Weak-coupling limit 

We follow Rcf. |lj and bosonize the Hamiltonian H. 
Since the bosonization of the XXZ chain is a standard 
procedure, we do not repeat the derivation of a bosonized 
Hamiltonian here. After performing the Jordan- Wigner 
transformation and taking continuum limit, we find that 
iJo reduces to a free-boson model, 



dx J 



(3) 



where n(a;) is a conjugate operator to the bosonic field 
(j){x): [(/){x),Il{y)] — iS{x — y). The spin wave velocity v 



is known to be w = (7r/20) sin 9, where = cos ^ A. The 
spins in the chain can be represented4ri terms of bosonic 
fields (t)ix) and 4>{x) (H = dij)/dx)\BM 
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Here x 



— ^ . ^^^.^ - - j and the c^'s are nu- 
merical constants. The lattice spacing is assumed to be 

I and (4b) is related 
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to A in the original Hamiltonian (R 



R 



Mil- 

— 1 cos 

27r V TT 

rib) 



1/2 



(5) 



With the Gaussian form of , we can immediately find 
the scaling dimensions of operators e"^"^ and e'"'^, both 
of which are a^/47r. Thus the dimensions of the stag- 
gered components of Sf and are (47ri?^)~^ and 7ri?^, 
respectively. 

From Eqs. ( |4a| ) and (^) the Kondo interaction term 
Hk becomes 



H 



(b) 
K 



imp 
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+ '^imp ( ^Fz + \Bz sm 
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where the couplings A's are proportional to Jk- Since 
the impurity spin is coupled to a single spin S^o in our 
model, we have backward Kondo scattering terms pro- 
portional to Xbz and A^j^. These terms do not appear in 
some models where 5'inip is coupled-ayjiametrically to two 
neighboring spins, say S'o and S'lEa^B These backscat- 
tering terms is an important ingredient of our model. The 
backward spinflip scattering term (oc \b±) has scaling 
dimension i:R^ and is always a relevant operator. This 
should be contrasted with the conventional Kondo prob- 
lem in 3D, where the Kondo interaction is a marginal op- 
erator of the form d(j)/dx. Therefore we conclude that the 
weak-coupling point {Jk = 0) is unstable for — 1 < A < 1 
independent of the sign of Jk, and the system always 
flows to a strong-coupling regime. This siijjation is quite 
similar to the Kondo effect in a TL liquid. □ To lowest or- 
der the scaling equation of the most divergent coupling 
Xbi_ is given by 



dXf 



dlogL 



= {l-TrR^)XB±, 



(7) 



where L is system size. We thus expect that the energy 
scale Tk at which the crossover from weak coupling to 
strong coupling occurs should be 



Tk oc \Xb±\ 



l/(l-7rfl'^) 



cx \Jk 



|l/(l-7rfl^) 



(8) 



for \Jk\ ^ J{— !)■ We-^dentify this energy scale with 
the Kondo temper ature.Ej 
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B. Strong-coupling limit for < A < 1 

Let us consider the strong-couphng limit where Jk 2> 
1. In this limit we first diagonalize Hk and treat the 
coupling between Sq and its neighbors (S^ii) as weak 
perturbations. The ground state of Hk is a spin singlet 
{So + S'imp = 0). In the limit Jk oo the system con- 
sists of the singlet and two decoupled semi-infinite chains 
(SICs). With very large but finite Jk, we derive effective 
interactions acting on the subspace of the singlet plus the 
SICs using \/ Jk expansionJH^ Second order perturbation 
yields 



1 



2Jk{1 + A) 
A2 



2J, 



K 



-SlSti + const. 



Higher order calculations also give the same form of in- 
teractions (and irrelevant operators). We now need to 
know the bosonization of these operators S±i at the 
boundaries of the SICs. This was discussed in detail 
by Eggert and AfHeckEj and we can simply borrow their 
results. The open-boundary condition implies that the 
phase field (t){x) is fixed to be some constant at a; = 0. 
To be specific, let us impose 0(0) = 0. The left- 
going field 4)l{x) = \/tt[4>{x) -\- 4>{x)] and the right-going 
field (pFiXx) — V^l'f'ix) — 0(a;)] are no longer indepen- 
dent. From these chiral fields we introduce two left-going 
fields: 



b<{x) = Q{-~x)(^l{x) - Q{x)(^r{-x), 



(10a) 
(10b) 



where Q{x) is a Heaviside step function. The field 4i>{x) 
defined on (—00,00) describes bosonic excitations in the 
SIC of the positive x region {Si: i > 0), and the other 
field 0<, also defined on (—00,00), describes excitations 
in the negative x region. Their commutation relations 
are [0> (x), 0> (y)] = [(?!)< (a;), 0< (y)] = -i7rsgn(x- y) and 
[(/)> (a;), 0< (y)] = 0. Their dynamics is governed by the 
Hamiltonian 



-ffsic = ^ I dx 



dx 



dx 



(11) 



With these fields the boundary spins can be written as 

d0>(O) 



S^ oc exp [±i2-\/7ri?0; 



Sti oc exp [±i2^R(t)< (0)] , S'Li 



oc 



dx ' 
rf 0<(O) 
dx 
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The scaling dimension of Sj^^ is 1 and that of is 

2ttR^. In general, the vertex operators e'""^ and e"*^ have 
dimension a^/2 at boundaries. We thus find that, among 
possible interactions generated by the 1/Jk expansions. 



SiS_i + Si S^i is most dangerous and has dimension 
AttR'^. This operator is irrelevant when < A < 1. 
Therefore we may conclude that, when the anisotropy 
parameter A of the host XXZ spin chain is < A < 1, 
the infrared stable fixed point corresponds to the limit 
Jk — 00, where the system is decoupled into a singlet 
and two semi-infinite XXZ spin chains; see Fig. 0. The 
singlet acts like an infinitely high potential barrier for ex- 
citations in the spin chain and effectively cuts it into two 
SICs. If the host spin chain is of finite length containing 
L spins and if we apply the periodic boundary condition, 
then its low-energy fixed point is an open spin chain con- 
sisting of L — 1 spins, in addition to a decoupled spin 
singlet formed from the impurity spin and a spin origi- 
nally in the host spin chainEj This strong-coupling fixed 
point is very similar to the one found for the Kondo effect 
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(9) in electronic TL liquids. I 
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FIG. 1. Schematic picture of renormalization to the 
strong-coupling fixed point where the XXZ chain is cut by 
the singlet. 

The above result is a natural generalization of the con- 
clusion of Eggert and Affieck to the case < A < 1. 
In their case the low-energy fixed point is a singlet plus 
decoupled two semi- infinite Heisenberg spin chains, and 
the leading irrelevant operator at the fixed point is a 
dimension 2 operator. Si ■ S-i. In our case the oper- 
ator Si SZi + Si S'^i has a smaller scaling dimension 
than S^S^i because of the absence of the SU(2) symme- 
try. Since its dimension AnR^ is in general noninteger, 
we may expect that it should give anomalous power-law 
temperature dependence to various quantities. 

The coupling to the Kondo impurity gives rise to an 
extra contribution to the specific heat and the spin sus- 
ceptibility, which we denote 5C and Sx- Their temper- 
ature dependence near the strong-coupling fixed point 
is determined by the leading irrelevant operator, Oi — 
gi {SiSZi + Si S^i^ , where gi is a coupling constant. 
To obtain leading temperaturojiependence we may use a 
perturbation expansion in Oi .E! 

We first estimate 6C. Up to second order, the change 
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in the free energy is given by 

SF^- dT(Oi(T)Oi(0)) 
Jo 

cx — / ar 



An ttTt J 



(13) 



where d = AttR^, (3 is inverse of the temperature T, 
Oi(t) = e^^^^^Oe^'^^^^'^ , and Tc is a cutoff to regularize 
the integral Note that there is no first-order contribu- 
tion of Oi to 5F. The low-temperature expansion of the 
integral in Eq. (|l^) for general d reads 



0/2 



ttTtc 



sinTrTr / 
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ld^i^Tr,)^'^-^Bill-d), 
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(14) 



where i3(a, 6) is the beta function. Note that any ir- 
relevant operator with dimension d > 3/2 generates a 
positive term. From these equations we get 



Ud- 



2d-2 



/(3-2d), 



<5Coc <^ riog(l/7rTT,), 
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d 
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in the low-temperature limit. Since d = AttR^, the 
boundary case d = | corresponds to A = l/\/2- When 

< A < l/\/2, 5C is proportional to T^''^'-^ with the 
exponent changing from to 1 as A varying from to 

1 / \/2- This anomalous power-law-behavior is reminiscent 
of the Kondo effect in TL liquids.Q The log correction ap- 
pears at A = l/\/2 when the dimension of the leading 
irrelevant operator becomes 3/2. This is mathematically 
the same as in the two-channel Kondo problem.Ej When 
l/\/2 < A < 1, the leading term of 5C is proportional to 
T. 

We next consider 5x- Here we need to distinguish two 
kinds of spin susceptibilities: one responding to a mag- 
netic field applied in the z direction and the other re- 
sponding to the one in the xy plane. We shall call them 
5xz and <5xt-i| respectively. Suppose we apply a magnetic 
field locallyEa only to 5'imp such that the perturbation, 



Hh — h,Si 



imp 



"'x'-'impi 



(16) 



is added to the Hamiltonian. Using the 1 / Jk expan- 
sion again, we can generate effective interactions in- 
duced by Hh in the Hilbert space of the singlet plus the 
SICs (Fig. |l|). From the symmetry we expect to have 
the following operators in addition to other less relevant 
ones: 0,a = h,{Sf + S^), 0^2 - hl{St SZ^ + S^S\), 
and Oh?, = hxiSf + S^i)- In terms of the bosonic 



fields they may be written as O/a oc hz[dx4>>{0) + 
5x</)<(0)], dh2 (X /i2cos{2V^i?[0>(O) - (^<(0)]}, and 
dh3 oc hx{cos[2y^R(j)>{0)] + cos[20ri?(/)<(O)]}, whose 
scaling dimensions are 1, AttR^, and 2'kR^. We can now 
estimate 5F induced by these operators using Eq. (|l^ ) 
and obtain 5xa — —d'^5F/dhWh=o- One point to be 
mentioned is that products of Oi and Oh2 can contribute 
a term h^T^^^ ~^ to SF, leading to a term proportional 
to T*^^ in Sxz- From these considerations, we con- 
clude that in the low-temperature limit Sx has the fol- 
lowing form: 



6xziT)-SxM^ 



6x±iT)^6x±{0)^T 




^ 0<A<l/x/2, 

A = l/V2, (17a) 
1/V2 < A < 1, 

< A < 1. (17b) 



We note that there is always a contribution propor- 
tional to coming from irrelevant operators. When 
< A < 1, the term T^^R^-^ might be difficult to ob- 
serve, because of its small coefficient (c>c AttR^ — 1), com- 
pared with the term. We also note that in general 
the zero-temperature limit of the susceptibility 6x{0) is 
of order 1/Tr-. 



C. Strong-coupling limit for 1 < A < 

When the parameter A in the host spin chain is in 
the range — 1 < A < 0, the dimension of the opera- 
tor S^SZi + Si is smaller than 1 and is relevant. 
This means that the open-boundary fixed point discussed 
in the previous subsection cannot be a low-energy fixed 
point when — 1 < A < 0. Both limits Jk and 
Jk oo in the original Hamiltonian Hq + Hk are un- 
stable. We thus need to find a nontrivial fixed point. 

Let us for the moment forget the singlet of S'imp and 
Sq, and concentrate on the rest of the spins. That is, we 
consider the two semi-infinite spin chains weakly coupled 
by a ferromagnetic exchange interaction H2'. 

oo 

H\ = ^ (SfSf_^i + SfSf^i + ASf Sf^i^j) 



1=1 

~X {SfSli + SlSli) , (18) 

where < A <C 1. We have dropped the irrelevant SfSZi 
term. Now we rotate Si (i > 0) around the z axis by tt, 
which changes the sign of A (—A A) in H\. We then 
apply RG transformation. Since SfS'Zi + SfS^^ is rele- 
vant, the coupling A grows as the energy scale decreases. 
The SfSti term is also generated in the course of the 
RG transformation. Thus, the two chains get coupled 
stronger at lower energy scale. We next consider the op- 
posite limit where the two chains are well connected but 
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one bond is slightly disturbed. This is described by the 
Hamiltonian, 



i=l 

oo 
i=l 

+ (1 - e^) {SfS^^, + S^Sl,) + (1 - e,.)ASlSl„ (19) 

where < £j_, £2 ^ 1. Bosonizing this Hamiltonian as in 
Sec. IIA, we find that the perturbations (oc e) give the 
spin-Peierlse operator sin [0(0) /i?] of dimension {4ttR'^)~^ 
and dimension 2 operators like {d(f>/dx)'^. Since they are 
irrelevant (ej_,£z — )■ in the low-energy limit), we re- 
cover a pure XXZ spin chain. It is tempting to assume 
that the RG trajectories starting from the unstable point 
describing two weakly coupled chains [Eq. (|l^)] contin- 
uously flow to the stable fixed point of the pure XXZ 
chain. Although we cannot prove it, we believe this is 
what actually happens. We note that this phenomenon is 
closely related to the well-known result that a backward- 
scattering potential is renormalized to zero fotfermions 
interacting with mutual attractive interactions. □ It is also 
similar to the "healing" of weak bonds which Eggert and 
Affleck found for the isotropic Heiaenberg chain with two 
symmetrically perturbed bonds .113 Coming back to the 
Hamiltonian H\, we conclude that its low-energy flxed 
point is a pure XXZ spin chain with the spins Si {i > 0) 
rotated around the z axis by tt. 
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FIG. 2. Schematic picture of renormalization to the sta- 
ble fixed point consisting of a singlet weakly coupled to a spin 
chain. The singlet looks "transparent" for low-energy excita- 
tions in the chain. 

We now return to our Kondo problem. What we have 
found so far is that (i) the Kondo coupling is a relevant 
operator at the weak-coupling point and leads to a sin- 
glet formation and that (ii) weakly coupled spin chains 
are renormalized to a strongly coupled single chain. Com- 
bining these two observations together, we propose the 
model schematically shown in Fig. |2| as a candidate for 



the low-energy fixed point. The model consists of the 
singlet of S^imp and Sq on top of the pure XXZ chain 
where spins are rotated as discussed in the last para- 
graph. An important point is that low-energy excitations 
are spin density fluctuations of long wave length in the 
chain and that for these low-energy excitations the sin- 
glet has essentially no effect. In other words, the singlet 
is "transparent" for them. At short-length scale there is 
a coupling between Sq and its neighbors {Si + S^i). 
We assume that, as far as low-energy physics is con- 
cerned, the singlet is rigid and can be broken only virtu- 
ally by the weak coupling of Sq to the spin chain. Thus, 
the stable fixed point may also be represented schemat- 
ically as in Fig. |3[ From the assumption of the rigid 
singlet, we can integrate out it to get effective interac- 
tions O2 = SfS'^i + SfSl^ and O3 = SfSi^ for the 
low-energy excitations in the spin chain. In the boson 
representation they are linear combinations of sin((/)/i?), 
{d4>/dxY, and (dcji/dt)'^, which are irrelevant operators 
for the spin chain with the parameter A in the range 
— 1 < A < 0. Hence the model is stable against weak 
perturbations, and we conjecture that the above model 
gives a correct picture of the strong-coupling fixed point 
for the case — 1 < A < 0. Although it is impossible to 
show analytically that the RG trajectories leaving from 
the unstable weak-coupling point reach this fixed point 
(Figs. ||and H), the numerical results we show in the next 
section provide good evidences for our picture. 



FIG. 3. Schematic picture of the low-energy fixed point. 

Assuming that our Kondo model is indeed renormal- 
ized to the strong-coupling fixed point of Fig. H, we can 
obtain leading temperature dependences of SC and 6x as 
in the last subsection. Since we know that a leading irrel- 
evant operator at the fixed point is among the operators 
sin[0(O)/i?], [d(j){0)/dx]^, and [d(j){0) / dt]^ , we find that 
the low-temperature behavior of 6C is given by Eq. ( p^ ) 
with d = (47ri?2)-i. We thus get 



T, 

SCoci Tlog(l/r), 

yl/(2irfl^)-2 



-1 < A < -1/2, 
A- -1/2, 
-1/2 < A < 0. 



(20) 



When a weak magnetic field is applied to S'imp, we 
obtain the operators Ohi — hz{Sl -t- Sti), Oh2 = 

hl{S+SZi + S^Sti), and 6,,3 = K{Sf + S-^i) after in- 
tegrating out the singlet. Since these operators are not 
boundary operators at the fixed point of our interest, the 
scaling dimensions of Oh2 and Ohi are different from the 
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open-boundary case. Here we use the bosonization for- 
mulas ( |4a| ) and (4b) and find that the dimensfons of Oh2 
and Oh3 are (47ri?^)~^ and ttR'^, respectively. We then 
obtain the following low-temperature behavior: 



oc <; r2iog(i/r), a 

^2|^^ R^T}2 



T'/il-6TTR'), 



i < A < 0, 

_i 

2 ' 



-1< A<-i, 



(21a) 
(21b) 



B. Numerical results for A — 0.5 

As a typical case of < A < 1 we have chosen A = 0.5. 
In this case R — 1/ a/Stt and v — 3^/3/4. With this choice 
we have computed lowest energy gap Eg for chains of 



L = 1 (mod 4 
shown in Fig. 



Numerical results of the finite-size gap is 
The energy gap is difference between the 



lowest energy in the sector S^^^ = and that in the sector 
Sfoj = 1. According to the RG analysis in Sec. IIB, the 
ground state of a sufficiently long chain is described as 
two decoupled chains, each having {L— l)/2 spins, plus a 
rigid spin singlet of Sq and Simp in between them. Note 
that (L — l)/2 = / is an even integer. 



D. Strong-coupling limit of the XY case (A = 0) 

We briefly comment on the low-energy fixed point for 
the XY case. Since this is exactly on the border of the 
two cases discussed in Sees. IIB and IIC, we naturally 
expect that a picture for the fixed point of the A = 
case should be something in between Figs. |l]and||. That 
is, the singlet of 5'imp and Sq does not completely cut 
the host XXZ spin chain into two pieces. The weakened 
connection between Si and S-i is not healed as in the 
negative A case. This is because at the open-boundary 
fixed point {Jk = oo) the operator S^SZi + S'^^ is a 
marginal operator. We expect that the impurity contri- 
bution to the specific heat and the susceptibilities have 
the following low-temperature limit: 

SC oc T, (22a) 

SxziT)-SxM'^T^, (22b) 
6x±iT)^log{l/T). (22c) 



III. RESULTS OF DMRG CALCULATIONS 

A. Numerical methods 

In this section we present our numerical results. The 
Hamiltonian we studied is Ho+Hk, Eqs. (0) and (|). The 
site index i in Eq. (|l|) runs from — Z to ^ — 1, and the total 
number of spins in the host XXZ chain is L = 21 + 1. We 
impose the open boundary condition at the left and right 
ends of the host XXZ chain. Using the DMRG method 
proposed by White,o we have calculated lowest energy 
gap and spin correlation functions in the ground state. 
In order to accelerate the numerical calculation, we haiie 
employed the improved algorithm proposed by White.Ej 
We have also used the finite system method to achieve 
high accuracy. Up to 100 states were kept for each block 
and the truncation error is typically 10~®. This error is 
directly related to the accuracy of energy. 
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FIG. 4. Energy gap Eg as a function of system size L for 
A — 0.5. The data points are the gap computed for L = 17, 
25, 33, 49, 65, 94, 129, 157, and 201 [L = 1 (mod 4)]. The 
dashed line represents the infinite-L limit, Egl — \^tv/2. 

To interpret finite-size scaling of the data, let us 
bosonize the two open XXZ chains of length I, follow- 
ing Refs. |l^ and |2^. The mode expansions of the phase 
fields are given by 

+ E (a«.e-*-* + at ^e'^^-*) , (23) 

~ ^ vt 

_^^^cosfc^ (24) 

where fc„ — T:n/l and the operators obey the commuta- 
tion relations [4>Qf,,Qu] = «<5^,iy and [am^,at^,] = 5m,„(5^,i. 
{fijV — I or r). The suffix / and r stand for the left {Sf. 
i < 0) and right {Si: i > 0) spin chain, respectively. The 
fields (t>i and (f>i {4>r and 0,.) are therefore defined in the 
negative (positive) x region, and 0/ -I- (pr and -I- 0^ cor- 
respond to and 4> in Eq. (^. Note that (/); and are 
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different fromji< and (/)>. Substituting Eqs. (^3|) and 
(^) into Eq. (y) yields the Hamiltonian of the fi chain 



\ n>0 / 



(25) 



Its energy eigenvalue and eigen functions are 



TTV 



n>0 



n 

Ji>0 



(26) 



-|0), (27) 



■ 



The constant Sf, 
of each 



where |0) is a vacuum (a„p|0) = 0) 
is nothing but a quantum number of total S 
chain. Since / is an even integer, can take integer val- 
ues only. Therefore, in the limit Jk — ^ oo, the ground 
state of the total system is the state with — TO„p — 
for n = I and r. The first excited states are fourfold 
degenerate and correspond to {Sf ,S^) — (±1, 0), (0, ±1) 



and 

by 



= 0. The energy gap in this limit is then given 



E, = ^27ri?2, 



(28) 



which equals ^Jii^j^l at A — 1/2. This gap value is 
shown as a dashed line in Fig. ^j. It is clear that all the 
curves in Fig. ^ are gradually approaching the dashed 
line as L increases. How the curves finally approach it in 
the L ^ oo limit is determined by the leading irrelevant 
operator Oi, whose explicit form we may take 



|i27ri? ^^(0,0)- ^,(0,0) |. (29) 



Oi oc cos ■ 



The correction to Eq. (|2^) due to the operator Oi cm, 
be obtained from lowest-order perturbation expansion.EII 
Since the degenerate first excited states \S[ = 1, — 0) 
and \Si = 0, = 1) have a nonzero matrix element. 



{S[ = 1, = 0\Oi\S[ = 0, = 1) 



oc (U| exp 
(X L-4-«' 



-27ri? 



1 



|0) 
(30) 



the degeneracy of these two states is lifted by an amount 
which scales as . The same is true for the other 

two degenerate states \S[ — ~1,S^ — 0) and \S[ — 
0,S^ = —1). On the other hand, the ground state en- 
ergy does not change in first-order perturbation. Hence 
we may expect that the leading correction to the energy 
gap should be proportional to ^-^Tri? ^ which goes to zero 
faster than the finite-size gap (oc L^^)- This L depen- 
dence is indeed observed in our numerical data shown 



in Fig. 1^. The data shows very clear power-law behav- 
ior with the exponent 4/3 — AttR^, in perfect agreement 
with the theory. This can be regarded as a numerical 
proof of the presence of the leading irrelevant operator 
with the scaling dimension AttR^ at the strong-coupling 
fixed point we discussed in Sec. IIB. We note that the 
energy gap Eg^'^ used in Fig. |5| is the one at Jj^ = oo, or 
equivalently, the finite-size gap of an XXZ spin chain con- 
taining I spins under the open boundary condition. The 
reason why we have used Eg^"^ rather than Eq. (^) is to 
reduce the effect of a bulk irrelevant operator, cos{2(j)/R), 
of dimension l/vri?^ = 3. 
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FIG. 5. Size dependence of the correction to the energy 
gap -Eg"', which is the gap calculated in the limit Jk = oo. 
The dashed line represents the theoretically predicted L"*^'^ 
dependence. 

Using the DMRG method, we have also calculated an 
equal-time two-point spin correlation function, (SfS^), 
in the ground state for L — 201 {S^^^ = 0). Accord- 
ing to our picture of the strong-coupling fixed point, the 
host XXZ spin chain is effectively cut by a singlet in the 
low-energy limit (Fig. 0). We naturally expect that cor- 
relations across the singlet should be much weaker than 
correlations within one of the decoupled chains. Our nu- 
merical results shown in Figs. ^ and ^ support this idea: 
A correlation function across the singlet show power-law 
dependence on i with an exponent larger than that for a 
pure XXZ chain {Jk = 0), 27ri?^ 

The exponents for (S^^Sf) and (S^^^^^^Sf) can be ob- 
tained from the following argument. First we consider 
{S'^^S~), which is equivalent to (S^^Sf). Since it van- 
ishes when the XXZ chain is completely decoupled, the 
nonzero contribution is due to the leading irrelevant oper- 
ator S^SZi + Si S'^-f^. To first order in Oi the correlator 
is 

(sl^s-) (X / dt{szdt)slM)i{stit)s-mr, (31) 
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where the averages ( ); and ( )r are evaluated for the 
ground state of each decoupled chain. Since the scaling 
dimension of the boundary operators S^i is 27r_R^ and 
that of S'^j is ttR^, we expect the correlator to scale as 

{Sl,Sn'^^-'^"^+\ (32) 

from which we get (S'^iSf) oc 1/i for A = 1/2. The 
results in Fig. ^ are consistent with this perturbative cal- 
culation. 




1 10 100 

i 



the whole system is a direct product of 15), which is the 
singlet wave function of S^imp and Sq, and the ground 
states of the left and right decoupled spin chains, which 
we denote |^) and \r). We calculate correlation function 
{S^pS~) to lowest order in the coupling between So and 
its neighboring spin, Sq S^: 

(SLpSr) ^ -L{S\S^^^\T){T\So){r\S+Sr\r) 

(x{~iyi-^^^\ (33) 

where \T) is a triplet state of Simp and So having exci- 
tation energy of order Jk- The exponent 37ri?^(= 1) is 
a sum of the dimensions of and S*". The data for 
Jk = 5 in Fig. |^ is in excellent agreement with the above 
calculation, although the data for Jk = 1 is curving, 
which we think is due to a crossover to the true scaling 
regime. 

C. Numerical results for A = —0.5 

Here we present the numerical results for negative A. 
Using the DMRG method, we have calculated finite-size 
gap and spin correlation functions for A = —0.5, where 
R = 1/VQtt and v = 3\/3/8. 



FIG. 6. Correlation between Sli and Sf calculated for 
A = 0.5 and L = 201. The dashed line corresponds to the 1/i 
decay obtained from the perturbative calculation. 
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100 



FIG. 7. Correlation between S^^p and Sf calculated for 
A = 0.5 and L — 201. The dashed line represents the ex- 
pected i~ behavior. 



The correlation between S'j^p and Sf 



can be calcu- 
lated using the 1/Jk expansion, which can be justified 
in the low-energy limit. At Jk = oo the ground state of 
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FIG. 8. Energy gap Eg as a function of system size L for 
A = -0.5. The data are taken for L = 17, 25, 33, 49, 65, 97, 
129, 201, and 301. The dashed line represents the infinite-L 
limit, EgL = VStv/S. 

Figure ^ shows the finite-size energy gap as a function 
of the system size for L — 1 (mod 4). As in the last 
section, the gap is defined as difference between the low- 
est energy in the sector S^^^. = and that in the sector 



1. We find that the normalized gap EgL increases 



for small Jk, while it decreases for large Jk- This is 
consistent with our picture of the renormalization flows 
(Fig. 1^). For small Jk the excitation gap is due to fluc- 
tuations of Simp weakly coupled to the host spin chain. 
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This coupling is renormalized and becomes stronger as 
we saw in Sec. IIA. For large Jk: on the other hand, the 
host XXZ chain is almost cut by a singlet, and the finite- 
size gap roughly corresponds to the singlet-to-triplet ex- 
citation energy in half chains. As L increases, or equiv- 
alently, as the energy scale decreases, the renormalized 
coupling between the almost decoupled chains becomes 
larger ("healing"), leading to the decrease of the normal- 
ized finite-size gap. It is clear that all the curves in Fig. ^ 
approach the dashed line EgL — ^/in/?) = 0.680 
which is the value one expects for a single XXZ chain 
of length L. Unlike in the case of A = 0.5, however, we 
have not been able to obtain information on the scaling 
dimension of a leading irrelevant operator from the nu- 
merical data. A log-log plot of jii^g — versus L did 
not give straight lines corresponding to power-law scal- 
ing. This would mean that the systems we have studied 
(L ^ 200) are not large enough. 
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FIG. 10. Correlation function {SUSt) for A = -0.5 and 
L = 201. The dashed line corresponds to the behavior. 
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We next discuss correlations between S'imp and Si. 
Since there is always a short-distance correlation be- 
tween Simp and Si + S-i, we expect (S'i^S'f ) oc ((Sf -I- 
S"i)S") with a smaller constant of proportion for larger 
Jk- Noting that Sf + S^i corresponds to the staggered 
component in a pure XXZ chain without the spin rota- 
tion of Si [i > 0), we conclude that (S'g„p5'f) cx i"^/^ 
and (Si^pSi) oc i""^ for large i. Our numerical results 
shown in Figs. 



: i ^ for large i. 
11 and O show exactly the feature dis- 



cussed above. Hence we conclude that the numerical re- 
sults support our picture of the low-energy fixed point. 



10" 



FIG. 9. Correlation function {SliSf) for A = -0.5 and 
L = 201. The dashed line corresponds to the decay. 

Next we show the results of correlation functions which 
we computed for the ground state of L = 201 system 
{S^^^ = 0) Figure | and |l^ show correlation functions 
of S^i and Si. The correlator (S't^Sf) is positive and 
decays like i"^/'^, while (S^^S^) is negative and decays 
much faster like These features are exactly what 

we expect from our picture of the low-energy fixed point 
(Figs. ^ and H). Since the spin chain is well connected, the 
correlation functions {S^.^S") should behave as in a pure 
XXZ chain without an impurity spin. That is, exponents 
of power-law decays should be the same as those in the 
pure chain, although amplitudes of the correlators will 
depend on Jk- From Eqs. (4a) and ( ^ ) we see that at 

- (-1)^6^2'^^'^, whose 
= 1/6. Hence (St.Sf) 
oc z^^/^, in agreement 



long distance Sf ~ d(j)/dx and 
scaling dimensions are 1 and ttR^ 
should decay as z""^ and LS^^Sf) 
with the numerical result .E3 
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FIG. 11. Correlation between Sf^p and Sf for A = -0.5 



and L 
havior. 



201. The dashed line corresponds to the i be- 
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FIG. 12. Correlation between Sf^^p and Sf for A = -0.5 
and L — 201. The dashed hne corresponds to the expected 
decay. 



in a pure spin chain {Jk = 0). These results are ex- 
plained successfully based on the RG analysis of the 
strong-coupling fixed point (Fig. |l|). For A = —0.5 we 
have found that the normalized gap approaches the value 
for a spin chain without the Kondo impurity. The corre- 
lation functions also show the same power-law behavior 
as in the pure spin chain. These results are consistent 
with our picture of the fixed point where the host spin 
chain remains as a single chain through healing of a cou- 
pling weakened by the singlet formation (Figs. || and|^). 
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IV. CONCLUSIONS 



In this paper we have studied the Kondo effect due to 
an extra spin coupled to a gapless XXZ spin chain. In 
our model the backward spinflip scattering is always a 
relevant perturbation. At low energy the impurity spin 
is screened by a spin in the host chain, and the char- 
acteristic energy scale, the Kondo temperature T^, has 
a power-law dependence on the Kondo coupling. From 
the perturbative RG analysis for various limits, we have 
deduced properties of strong-coupling, low-energy fixed 
points. In the antiferromagnetic side (0 < A < 1) 
the host XXZ chain is cut by the singlet into two sep- 
arate chains. On the other hand, in the ferromagnetic 
side (— 1 < A < 0) the singlet does not harm the host 
spin chain in the low-energy limit. This may be under- 
stood qualitatively by mapping the problem to a spinless 
fermionic system using the Jordan- Wigner transforma- 
tion. The fermions have mutual repulsive (attractive) 
interactions in the antiferromagnetic (ferromagnetic) re- 
gion. The singlet may then be viewed as an impurity 
potential for the fermions, which can be a relevant or ir- 
relevant perturbation, depending on the sign of the mu- 
tual interactions. Employing the known result for the 
spinless fermion system,Ll we can argue that the host spin 
chain is cut into two pieces in the antiferromagnetic case 
whereas in the other case the singlet does not affect the 
low-energy properties of the spin chain. 

We have used the powerful DMRG method to nu- 
merically compute finite-size energy gaps and correla- 
tion functions. The numerical results are consistent with 
the RG analysis. For A — 0.5, the normalized gap ap- 
proaches the value for a chain of half length, and the 
correlation function across Sq decays much faster than 
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